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ABSTRACT 

It has been widely shown that the cosmological parameters and dark energy 
can be constrained by using data from type-la supernovae (SNe la), the cosmic 
microwave background (CMB) anisotropy, the baryon acoustic oscillation (BAO) 
peak from Sloan Digital Sky Survey (SDSS), the X-ray gas mass fraction in 
clusters, and the linear growth rate of perturbations at 2; = 0.15 as obtained 
from the 2dF Galaxy Redshift Survey. Recently, gamma-ray bursts (GRBs) 
have also been argued to be promising standard candles for cosmography. In this 
paper, we present constraints on the cosmological parameters and dark energy by 
combining a recent GRB sample including 69 events with the other cosmological 
probes. First, we find that for the ACDM cosmology this combination makes the 
constraints stringent and the best fit is close to the fiat universe. Second, we fit 
the fiat Cardassian expansion model and find that this model is consistent with 
the ACDM cosmology. Third, we present constraints on several two-parameter 
dark energy models and find that these models are also consistent with the ACDM 
cosmology. Finally, we reconstruct the dark energy equation-of-state parameter 
w{z) and the deceleration parameter q{z). We see that the acceleration could 
have started at a redshift from Zt = OAO^q'^ to zt — O.GS^qoj. This difference 
in the transition redshift is due to different dark energy models that we adopt. 
The most stringent constraint on w{z) lies in the redshift range ^ ~ 0.3 — 0.6. 
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1. Introduction 

The traditional cosmology has been revolutionized by modern observational techniques 
in distant Type la supernovae (SNe la) (Riess et al. 1998; Perlmutter et al. 1999), cosmic 
microwave background (CMB) fluctuations (Bennett et al. 2003; Spergel et al. 2003, 2006), 
and large-scale structure (LSS) (Tegmark et al. 2006). These observations suggest that the 
composition of the universe may consist of an extra component such as dark energy or the 
equations governing gravity may need a variation to explain the acceleration of the universe 
at the present epoch. 

SNe la have been considered as astronomical standard candles and used to measure 
the geometry and dynamics of the universe. However, since it is difficult to observe SNe la 
at redshift z > 1.7, this measurement has been carried out only for the z < 1.7 universe. 
Recently, it was shown that GRBs may be complementary to the SN cosmology for three 
reasons. First, GRBs are the most powerful explosive events at cosmological distances and 
in particular long-diiration GRBs originate from the core collapse of massive stars. So 
GRBs would be detectable out to very high redshifts when the core collapse of the first 
stars occur (Ciardi & Loeb 2000; Lamb & Reichart 2000; Bromm & Loeb 2002, 2006). In 
fact, the farthest burst detected so far is GRB 050904, which is at 2; = 6.295 (Kawai et 
al. 2006). Thus, GRBs could provide a much longer arm for measuring changes in the 
slope of the Hubble diagram than do SNe la. Second, gamma-ray photons suffer from no 
dust extinction when they propagate to us, so the observed gamma-ray flux is a direct 
measurement of the prompt emission energy. Third, there have been extensive discussions 
on relations between the spectral and temporal properties and some of these relations have 
been shown to be promising standard candles for cosmography. Schaefer (2003) derived the 
luminosity distances of 9 GRBs with known redshifts by using two quantities (the spectral 
lag and the variability) as luminosity calibrators and gave a constraint on the mass density 
JIm- Ghirlanda et al. (2004a) found a tight relation between coUimation-corrected energy 
and the local-observer peak energy (i.e., the so-called Ghirlanda relation). This relation 
may be physically understood as due to the viewing angle effect of an annular jet (Levinson 
& Eichler 2005) or Comptonization of the thermal radiation flux that is advected from the 
base of an outflow (Rees & & Mcszaros 2005; Thompson et al. 2006). Assuming that some 
physical explanation (e.g., the understandings mentioned above) comes into existence, Dai, 
Liang & Xu (2004) used the Ghirlanda relation to constrain the cosmological parameters 
and dark energy. Since then, a lot of work in this so-called GRB cosmology fleld has been 
pubhshed (Ghirlanda et al. 2004b; Di Girolamo et al. 2005; Firmani et al. 2005; Friedman 
& Bloom 2005; Lamb et al. 2005; Liang & Zhang 2005, 2006; Xu, Dai & Liang 2005; Wang 
& Dai 2006a; Li et al. 2006; Su et al. 2006; Schaefer 2007; Wright 2007). Very recently, 
Schaefer (2007) used 69 GRBs and five relations to build the Hubble diagram out to 2; = 6.60 
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and discussed the properties of dark energy in several dark energy models. He found that 
the GRB Hubble diagram is consistent with the concordance cosmology. Besides SNc la 
and GRBs, the other observations such as the shift parameter of CMB (Spergel et al. 2003, 
2006), the baryon acoustic peak from Sloan Digital Sky Survey (SDSS) (Eisenstein et al. 
2005), the X-ray gas mass fraction in clusters (Allen et al. 2004), the perturbation growth 
rate from 2dF Galaxy Redshift Survey (Hawkins et al. 2003), and the weak lensing (e.g., 
Schimd et al. 2007) have been used to constrain cosmological parameters and explore the 
properties of dark energy. 

It is of growing interest that dark energy is reconstructed in a model-independent way 
to investigate the evolution of the deceleration parameter q[z) and the dark-energy equation- 
of-state parameter w{z) (Alam et al. 2004; Vircy et al. 2005; Gong & Wang 2007; Alam et 
al; 2007). Evolving dark energy models had been shown to satisfy the data from SNe la. To 
reconstruct q{z) and w{z), Gong & Wang (2007) used the new "Gold" sample of SNe and 
data of SDSS and CMB, while Alam ct al. (2007) adopted the new "Gold" SN sample, the 
SNLS sample, and data of SDSS and CMB. It is found that the result is strongly dependent 
on the matter density VLm- The transition redshift Zt ~ 0.2 was found in reconstruction of 
q{z) (Virey et al. 2005; Shapiro & Turner 2006; Gong & Wang 2006). Previous investigations 
in the construction of w{z) show that the stringent constraint on w{z) is in the redshift range 
^ ~ 0.2 - 0.5 (Alam et al. 2004; Gong & Zhang 2005). 

In this paper we use GRBs and the other observational data to measure the cosmological 
parameters and the nature of dark energy. We also reconstruct q{z) and w{z) out to 2; > 
6.0 using these observational datasets, explore the transition redshift and constrain w{z). 
Recently, Su et al. (2006), Li et al. (2006), and Wright (2007) combined GRBs with some 
other cosmological probes to constrain the ACDM cosmology, the constant w model, and 
the dark energy model of w{z) — WQ-\-Wa{^ — a) (where a is the scale factor of the universe), 
respectively. In their papers, these authors adopted the distance modulus and its error (of 
a GRB) calculated for the concordance cosmology or the dynamical dark energy model of 
w{z) = —1.31 + lA8z, which were presented by Schacfcr (2007). In addition, Li et al. (2006) 
used the Markov Chain Monte Carlo technique to carry out global fitting. Here we use the 
observational data (e.g., time lag, variability, spectral peak energy -Epeak, minimum rise time) 
of GRBs to make a simultaneous fit of five correlations in any given cosmology, and consider 
more other cosmological probes and more dark energy models. The structure of this paper is 
arranged as follows: in section 2, we introduce GRBs and the other cosmological probes and 
describe our analytical methods. The constraints on the cosmological parameters and dark 
energy are presented in section 3. In section 4, we reconstruct w{z) and q{z). In sections 5, 
we summarize our findings and present a brief discussion. 
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2. Observational data and Analysis Methods 



2.1. Type la Supernovae (SNe la) 

Riess et al. (2004) reanalyzed the SN la dataset. They considered 14 new high-redshift 
events observed by the Hubble Space Telescope (HST). This led to a sample known as the 
"Gold" sample containing 157 SNe la. Recently Riess et al. (2007) added 25 SNe la to this 
sample. The final sample now consists of 182 SNc la. The observations of SNe la provide the 
currently most direct way of probing the dark energy at low-to-medium redshifts because 
the used luminosity distance is directly related to the expansion history of the universe, that 
is, 

cH^\l + z)(-l]fc)-^/2gin[(-(]fc)V2/] < 0, 

dL^ < cHo\l + z)I Qjk = 0, (1) 



1 — f^M — ^ 



I'DE, 



where 

(2) 
(3) 

(4) 

where w{z) is the equation-of-state parameter for dark energy and dL is the luminosity 
distance. With dL in units of megaparsecs, the predicted distance modulus is 



7= / dz/E{z), 
Jo 

E{z) = [(1 + zfnM + f{z)nDE + (1 + zY^kV^^ 

f{z) = exp ' 



= 51og(dL) + 25. 



(6) 



The likelihood functions for the parameters VLm and VLde can be determined from statis 
tics, 

X [rlo, i^M, "deJ = > ^ o , — o ' 

' * -4- 



(7) 



where a^, is the dispersion in the supernova redshift (transformed to distance modulus) due 
to a peculiar velocity, /i.o,i is the observed distance modulus, and cTjuq . is the uncertainty 
in the individual distance modulus. The confidence regions in the — ^de plane can be 
found through marginalizing the likelihood functions over Hq (i.e., integrating the probability 
density p oc exp(— x2/2) for all values of Hq). 
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2.2. Gamma-Ray Bursts (GRBs) 



GRBs can be detected out to very high redshifts (Ciardi & Loeb 2000; Lamb & Reichart 
2000; Bromm & Loeb 2002, 2006). They can bridge up the gap between the nearby SNe la 
and the distant CMB anisotropy. Schaefer (2007) comphed 69 GRBs to make simultaneous 
uses of five luminosity indicators, which are relations of riag — I/, V — L, Epeak — L, £'peak — E^, 
and trt — L. Here the time lag (riag) is the time shift between the hard and soft light curves, 
L is the luminosity of a GRB, the variabihty F of a burst denotes whether its hght curve 
is spiky or smooth and V can be obtained by calculating the normalized variance of an 
observed light curve around a smoothed version of that light curve (Fenimore & Ramirez- 
Ruiz 2000), Epeak is the peak energy in the z/Fj, spectrum, E^ = {1 — cos9j)Eiso is the 
collimation-corrected energy of a GRB, and the minimum rise time (trt) in the gamma-ray 
light curve is the shortest time over which the light curve rises by half of the peak flux of 
the pulse. We make a simultaneous fit to these five relations for any fixed cosmology. We 
perform a linear regression analysis to find a relation between observational quantities. After 
obtaining the distance modulus of each burst using one of these relations, we use the same 
method as Schaefer (2007) to calculate the real distance modulus. 



where the summation runs from 1 — 5 over the relations with available data, //.j is the best 
estimated distance modulus from the i-th relation, and cr^. is the corresponding uncertainty. 
The uncertainty of the distance modulus for each burst is 



Fig.l shows the Hubble diagram from the new "Gold" SNIa sample and 69 GRBs. The 
combined Hubble diagram is consistent with the concordance cosmology. GRBs can build 
the Hubble diagram out to 2; > 6.0 (Schaefer 2007). The GRB Hubble diagram is well- 
behaved and describes the shape of the Hubble diagram at high redshifts. When calculating 
constraints on cosmological parameters and dark energy, we do not care about the slopes 
of the five relations because we have marginalized these parameters (Schaefer 2007). The 
marginalization method is to integrate over some parameter for all of its possible values. We 
also marginalize the nuisance parameter Hq. The value is 




(8) 




(9) 




(10) 



where //fit,i ^-nd (7^^^ . are the fitted distance modulus and its error. 
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2.3. Cosmic Microwave Background (CMB) 

Observations of the CMB anisotropy provide us with very accurate measurements, which 
may be used to gain insight about dark energy and cosmological parameters (Spergel et al. 
2006). We may make use of the 3-year WMAP results to get the shift parameter (Wang & 
Mukherjee 2006) 

where E{z) = H{z)/Hq and the function sinn(x) is defined as sinn(x) = sin(2;) for a closed 
universe, sinn(x) = sinh(x) for an open universe and sinn(x) = x for a flat universe. To 
calculate the last scattering redshift z^^, we adopt Jl^^^ — 0.024 and Qj^h'^ — 0.14 ± 0.02. To 
calculate zis, we consider a fitting function: 

zi, = 1048[1 + 0.00124(O6/i2)-0.738][;L ^ g.i^Mh^n (12) 

where the quantities gi and g2 are defined as gi = 0.078{^bh^y^-'^^^[l + 39.5{ilbh^)^''^^^]~'^ 
and g2 = 0.56[1 + 21.1{nbh^y-^^]-^ respectively (Hu & Sugiyama 1996). The value is 

XCMB — g ■ ^^"^^ 



2.4. Baryon Acoustic Peak from SDSS 

It is well known that the acoustic peaks in the CMB anisotropy power spectrum can 
be used to determine the properties of perturbations and to constrain cosmological param- 
eters and dark energy (Spergel et al. 2003). The acoustic peaks occur because the cosmic 
perturbations excite sound waves in the relativistic plasma of the early universe (Peebles 
& Yu 1970; Holtzmann 1989). Because the universe has a fraction of baryons, the acoustic 
oscillations in the relativistic plasma would be imprinted onto the late-time power spectrum 
of the nonrelativistic matter (Peebles & Yu 1970; Eisenstein & Hu 1998). The acoustic 
signatures in the large-scale clustering of galaxies can also be used to constrain cosmological 
parameters and dark energy by detection of a peak in the correlation function of luminous 
red galaxies in the SDSS (Eisenstein et al. 2005). This peak can provide a "standard ruler" 
with which the cosmological parameters and dark energy are measured. We use the value 



rsmn 



E{z,)\nk\ V Jo E{z) 



1/3 



(14) 
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measured from the SDSS data to be A = 0.469(0.95/0.98)-°-^^ ± 0.017, where Zi = 0.35. 
The value is 

_ {A- 0A69f 
Xbao - 0.0172 ■ ^ ^ 



2.5. X-ray Gas Mass Fraction in Clusters 

Since clusters of galaxies are the largest virialized systems in the universe, their matter 
content is thought to provide a sample of the matter content of the universe. A comparison 
of the gas mass fraction, /gas = Mgas/Mtot, as inferred from X-ray observations of clusters 
of galaxies to the cosmic baryon fraction can provide a direct constraint on the density 
parameter of the universe Qm (White et al. 1993). Moreover, assuming the gas mass 
fraction is constant in cosmic time, Sasaki (1996) showed that the /gas measurements of 
clusters of galaxies at different redshifts also provide an efficient way to constrain other 
cosmological parameters decribing the geometry of the universe. This is based on the fact 
that the measured /gas values for each cluster of galaxies depend on the assumed angular 
diameter distances to the sources as /gas oc (D^)'^/^ xhe true, underlying cosmology should 
be the one which makes these measured /gas values invariant with redshift (Sasaki 1996; 
Allen at al. 2004). Using the Chandra observational data, Allen et al. (2004) have got the 
/gas profiles for the 26 relaxed clusters. These authors used the 26-cluster data to constrain 
cosmological parameters. They found = 0.245^^:^3? and Qa = OMtlH in the ACDM 
cosmology. This database has also been used to constrain the generalized Chaplygin gas 
model (Zhu 2004) and the braneworld cosmology (Zhu and Alcaniz 2005). We will combine 
this probe in our analysis. Following Allen et al. (2004), we calculate the value as 

/ 26 r /■SCDM/„\ _ f l2 
2 _ I L-'gas V^V ./gas, ij 

\ 1=1 /gas, 1 



/Qb/i^- 0.0233 V //i- 0.72 V /6- 0.824V , ^ 

+ [ 0.0008 ) +l^:o^J ^[-^mrj ' ^^'^ 

where f^^^^{z) = bQb/[{'^ + 0.19y/h)QM] x [df^^{z)/dT'^{z)Y-^, /gas,i is the observational 
baryon gas mass fraction and fe is a bias factor motivated by gas-dynamical simulations which 
suggest that the baryon fraction in clusters is slightly lower than for the universe as a whole. 
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2.6. Perturbation Growth Rate from 2dF Galetxy Redshift Survey 

The clustering of galaxies is determined by the initial mass fluctuations and their evolu- 
tion. We can set constraints on the initial mass fluctuations and their evolution by measuring 
the galactic two-point correlation function. The 2dF galaxy redshift survey measured the 
two point correlation function at the redsift oi z — 0.15. Hawkins et al. (2003) measured the 
redshift distortion parameter P — 0.49 ± 0.09. This result can be combined with the linear 
bias parameter b — 1.04±0.11. So the growth factor g a,t z — 0.15 is g — bx P — 0.51 ±0.11. 
Theoretically, this growth factor is cosmology-dependent. Thus, the measurement of the 
perturbation growth rate (PGR) g{z — 0.15) can be used to calculate x^'- 



_ {g-0.5ir 

which constrains the cosmological parameters and dark energy. 



3. Constraints on Cosmological Parameters and Dark Energy 

Using the datasets of the above observational techniques, we measure cosmological pa- 
rameters and dark energy. We can combine these probes by multiplying the likelihood 
functions. The total value is 

Xtotal ~ XSN + XgRB + XCMB + XbAO + Xgas ' XPGR \^°) 



3.1. The ACDM Cosmology 

The luminosity distance in a Friedmann- Robertson- Walker (FRW) cosmology with mass 
density Qm and vacuum energy density (i.e., the cosmological constant) fl\ is (Carroll, Press 
& Turner 1992) 

dL = c{l + z)H^^\nk\~^/^smn{\nk\^/^ 

X [ dz[{l + zf{l + nMz)-z{2 + z)nA\-^^''}. (19) 
Jo 

We use the datasets discussed above to constrain cosmological parameters. Fig.2 shows the 
la contours plotting in the Qm — ^a plane. The thick black line contour from all the datasets 
shows = 0.27 ± 0.02 and Qa = 0.73 ± 0.08 (la) with Xmin = 270.60. The red contour 
shows a constraint from 69 GRBs, and for a flat universe, we measure 0,m — 0.341q;°q (1<t), 
which is consistent with Schaefer (2007). Because the thin sohd line in Fig.2 represents a 
flat universe, our result from all the datasets favors a flat universe. 
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3.2. The Cairdassian Expansion Model 

The Cardassian expansion models (Freese & Lewis 2002) involve a modification of the 
Friedmann equation, which allows an acceleration in a fiat, matter-dominated cosmology. 
We assume that the Cardassian expansion model is (Freese & Lewis 2002; Zhu et al. 2004) 

H'-^-^ip + Cp^. (20) 

This modification may arise from embedding our observable universe as a (3+l)-dimensional 
brane in extra dimensions or the self-interaction of dark matter. The luminosity distance in 
this model is 

dL = cHQ\l + z) f dz[{l + zf^M + {I - ^m){1 + zf'']-^/^ (21) 

Fig.3 shows constraints on Qm and n. The sohd contours are obtained from all the datasets. 
Prom this figure, we have Qm — 0.28 ±0.02 and n — 0.02lQ;Qg at the la confidence level with 
Xmin — 272.52. This result is consistent with the ACDM cosmology. 



3.3. The w{z) = wq Model 

We consider an equation of state for dark energy 

w{z) = Wq. (22) 
In this dark energy model, the luminosity distance for a fiat universe is (Riess et al. 2004) 

dL = cHQ\l + z) f dz[{l + zfnM + {l-^M){l + zf^^+'"'Y^/^ (23) 
Jo 

Fig.4 shows the constraints on wq versus Jl^ in this dark energy model from all the datasets. 
Prom this figure, we have Qm = 0.31 ± 0.03 and wo = -O.QSIq^s (Ifj) with Xmin = 272.23. 



3.4. Two-Pctrameter Dark Energy Models 

Using the parameterization 

/ X wiz , 

w{z) = wo + -—, (24) 
1 + z 

the luminosity distance is calculated by (Chevallier & Polarski 2001; Linder 2003) 

dL = cHq\1 + Z) r dz[{l + z)^nM + (1 - ^m){1 + ^)3(l+"'0+«'i)g-3«;i./(l+.)]-l/2_ ^^5) 

Jo 
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Fig.5 shows the constraints on wq versus Wi in this dark energy model. The sohd contours are 
obtained from all the datasets and we find xLin — 273.25. Wq = — LOSlggg and Wi = 0.84lQg2 
(Icr) for the prior of = 0.30. We also assume this prior in the following analysis. 

Jassal, Bagla and Padmanabhan (2004) modified the above parameterization as 

w{z) = wo+ (i^^^)2 - (26) 

This equation can model a dark energy component which has a similar value at lower and 
higher redshifts. The luminosity distance is 

dL = CH^\1 + Z) r dz[{l + zf^M + (1 - Om)(1 + ^)3(l+-o)g3^i.V2(l+.)^]-l/2_ ^3?) 
^0 

Constraints on wq and w\ are presented in Fig.6. Prom this figure, we find Xmin — 272.07, 
wq — — 1.36lo;48 and wi ~ 3.32^182 (la) from all the datasets (blue contour). 

The third dark energy model that we consider is (Alam et al. 2003) 

1 + ^^ + 2^2(1 + ^) , . 

"^^^ = ^ ^) ^''^ 

where ^de{z) is defined as 

^de{z) = Ai{l + z) + A-2{1 + 2)2 + 1 - Qm - A - A2. (29) 

Fig. 7 shows the constraints on Ai versus A2 in this dark energy model. The solid contours are 
obtained from all the datasets and we find xLn = 273.95, Ai = -QAStiH and A2 = 0.221^3^ 
(la). 



4. Reconstruction of w{z) and q{z) 

Many dark energy models have been proposed (Copcland et al. 2006; Bludman 2006 
for a recent review) and we have fitted these models using the observational data in the last 
section. We now explore the properties of dark energy in a model-independent way (Sahni et 
al. 2006 for a review). In the following we reconstruct dark energy to find new information 
about dark energy from most of the recent datasets. The method to reconstruct directly 
properties of dark energy from observations in a quasi-model independent method has been 
discussed (Alam et al. 2004; Gong & Wang. 2007; Alam et al. 2007). We determine the 
dark energy equation of state based on 

2/1 I ~^ dlnH _ -1 
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The deceleration parameter 



q{z)^{l + z)H 



- 1. 



(31) 



dz 



We consider the first ansatz 



H{z) = //o[(l + zf^M + (1 - ^m){1 + ^)3(l+-o+«.i)g-3«,i./(l+.)]l/2]^ 



(32) 



which is in fact equivalent to the parameterization equation (24). The evolution of w{z) is 
plotted in Fig. 8. It is easy to see that the errors of the constraint on the equation of state 
become larger with redshift. The stringent constraint on w{z) happens z — 0.3 ~ 0.7. 
Using the GRB data, we can reconstruct w{z) out to z ~ 6.0 in the bottom panel. The 
evolution of q{z) is plotted in Fig. 9. We can see that the transition redshift at which the 
expansion of the universe was from deceleration {q{z) > 0) to acceleration {q{z) < 0) is 
zt = 0.57^0 07 (Ic'')- This result is consistent with Riess et al. (2004) and Wang & Dai 



which is in fact equivalent to the parameterization equation (26). The evolution of w{z) is 
plotted in Fig.lO. The stringent constraint on w{z) happens at 2; = 0.2 ~ 0.35. Using the 
GRB data, we can reconstruct w{z) out to 2; ~ 6.0 in the bottom panel. We find that the 
constraint on w{z) is also stringent around z = 4.0 ~ 5.0. The evolution of q{z) is plotted 
in Fig. 11. We can see that the transition redshift is zt — 0.40l2!o8 (l*^)- 

We consider the third ansatz 



H{z) = Ho[{l + zf^M + Ai(l + z) + A2{1 + zf + l-nM- A^- A^Y'^ (34) 



which is in fact equivalent to the parameterization equation (28). The evolution of w{z) is 
plotted in Fig. 12. The stringent constraint on w{z) happens at 2; = 0.35 ~ 0.55. Using the 
GRB data, we can reconstruct w{z) out to 2; ~ 6.0 in the bottom panel. We find that the 
constraint on w{z) becomes stringent around z ~ 6.0. The evolution of q{z) is plotted in 
Fig. 13. ^From this figure, we can see that the transition redshift is zt = 0.65lo;o5 (1"^)- 



(2006a, 2006b). 



We consider the second ansatz 



H{z) = H^[{1 + zfVLM + (1 - Om)(1 + ^)3(i+-o)e3«„.V2(i+.)Y/2^ 



(33) 



5. 



Conclusions and Discussion 



In this paper, we have presented the constraints on the cosmological parameters and dark 
energy by combining a recent GRB sample including 69 events with the 182 SNe la, CMB, 
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BAO, the X-ray gas mass fraction in clusters and the hnear growth rate of perturbations at 
z — 0.15 as obtained from the 2dF galaxy redshift survey. Wc found that the mass density 
of the universe is = 0.27 ± 0.02 and Oa = 0.73 ± 0.08 (Ifj) in the ACDM cosmology 
This result is well consistent with a flat universe. We also found that = 0.28 ± 0.02 and 
n = 0.02^Q Qg (la) in the flat Cardassian expansion model. We fltted several dark energy 
models. Finally, we reconstructed the dark energy equation-of-state parameter w{z) and the 
deceleration parameter q{z). We found that the the cosmic acceleration could have started 
between the redshift zt — 0.40lo;o8 — O.GSlo^os (l"^)- The stringent constraints on 

w{z) lie in the redshift range 2; ~ 0.3 — 0.6. 

Based on our analysis, it can be seen that the preferred cosmological model is the 
flat ACDM cosmology because of a small minimum value, Xmm — 270.60. The other 
models such as the Cardassian expansion model, the flat constant w model, and three two- 
parameter dark energy models can also flt all the datasets because the minimum values 
in these models vary only from Xmin — 272.23 to Xmin — 273.95. Thus, we cannot reject any 
of these models. 

It is well known that the cosmological constant suffers from the "flne tuning" problem 
and the coincidence problem (Zeldovich 1968; Weinberg 1989). In this paper, therefore, 
we have considered alternative possibilities, e.g., the Cardassian expansion model, the flat 
constant w model, and three two-parameter dark energy models. As we have shown, all 
the alternative models can be reduced to the flat ACDM cosmology at the la confidence 
level. So one needs more new observed data to distinguish between these models. New 
observations would be expected to improve the current constraints and test the flat ACDM 
model. GRBs appear to be natural events to study the universe at very high redshifts. The 
forthcoming GLAST will accumulate more GRB data, and in particular, its combination 
with Swift would lead to stronger constraints on high-redshift properties of dark energy. 
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Fig. 1. — Hubble diagram of new 182 SNe la (filled circles) and 69 GRBs (open circles). The 
solid line is calculated for a flat cosmology: = 0.27 and Qa = 0.73. 
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Fig. 2. — The la joint confidence contours for {Qm, ^a) from tfie observational datasets. Tlie 
tliick black line contour corresponds to all the datasets. The blue contour corresponds to 26 
galaxy clusters. The red contour corresponds to 69 GRBs. The yellow contour corresponds to 
the CMB shift parameter. The violet contour corresponds to 182 SNe la. The orange contour 
corresponds to BAO. The purple contour corresponds to 2dF Galaxy Redshift Survey. The 
thin solid line represents a flat universe 
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Fig. 3. — The la and 2a joint confidence contours for {flM,n) from all the observational 
data. 



- 19 - 




Fig. 4. — The la, 2a and 3a joint confidence contours for {Qm, wq) from all the observational 
datasets. 
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Fig. 5. — The la and 2a joint confidence contours of from all the observational data in the 
w{z) — wq + wiz/ (1 + z) model. 
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Fig. 6. — The Icr joint confidence contours from tlie observational datasets in tlie w{z) = 
wo + wiz/{l + z)'^ model. The blue contour corresponds to all the datasets. The cyan contour 
corresponds to 69 GRBs. The black contour corresponds to the CMB shift parameter. 
The green contour corresponds to 182 SNe la. The dashed contour corresponds to BAO. 
The yellow contour corresponds to the perturbation growth rate from 2dF Galaxy Redshift 
Survey. 
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Fig. 7. — The la and 2a joint confidence contours for (741,742) from all the observational 
data in the model of equation (28). 
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Fig. 8. — The evolution of w{z) by fitting the model w{z) = Wq + Wiz/{1 + z) to all the 
observational data. The solid line represents the reconstructed w{z). The shaded region 
shows the la error. We can constrain the evolution of w{z) up to z > 6.0 using GRBs 
(bottom panel). 
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Fig. 9. — The evolution of q[z) by fitting the model w{z) = Wq + Wiz/{1 + z) to all the 
observational data. The solid line is plotted by using the best fitting parameters. The 
shaded region shows the la error. We can reconstruct q{z) up to z > 6.0 using GRBs 
(bottom panel). 
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Fig. 10. — Same as Fig.8 but fitting the model w{z) = Wo + Wiz/{1 + z)"^ to all the observa- 
tional data. 
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Fig. 11. — Same as Fig.9 but fitting the model w{z) = Wo + Wiz/{1 + z)"^ to all the observa- 
tional data. 
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Fig. 12. — Same as Fig.8 but fitting the model of equation (28) to all the observational data. 
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Fig. 13. — Same as Fig.9 but fitting the model of equation (28) to all the observational data. 



